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Extrinsic point defects in aluminum antimonide 
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We investigate thermodynamic and electronic properties of group IV (C, Si, Ge, Sn) and group VI (O, S, 
Se, Te) impurities as well as P and H in aluminum antimonide (AlSb) using first-principles calculations. To 
this end, we compute the formation energies of a broad range of possible defect configurations including defect 
complexes with the most important intrinsic defects. We also obtain relative scattering cross strengths for 
these defects to determine their impact on charge carrier mobility. Furthermore, we employ a self-consistent 
charge equilibration scheme to determine the net charge carrier concentrations for different temperatures and 
impurity concentrations. Thereby, we are able to study the effect of impurities incorporated during growth and 
identify optimal processing conditions for achieving compensated material. The key findings are summarized 
as follows. Among the group IV elements, C, Si, and Ge substitute for Sb and act as shallow acceptors, while 
Sn can substitute for either Sb or Al and displays amphoteric character. Among the group VI elements, S, Se, 
and Te substitute for Sb and act as deep donors. In contrast, O is most likely to be incorporated as an interstitial 
and predominantly acts as an acceptor. As a group V element, P substitutes for Sb and is electrically inactive. C 
and O are the most detrimental impurities to carrier transport, while Sn, Se, and Te have a modest to low impact. 
Therefore, Te can be used to compensate C and O impurities, which are unintentionally incorporated during the 
growth process, with minimal effect on the carrier mobilities. 

PACS numbers: 61.72J-, 72.10.Fk, 71.55.Eq, 61.72.U-, 71.15.Mb 



I. INTRODUCTION 

Applications in nuclear and medical imaging as well as 
homeland security have recently led to a resurgence of re- 
search in the field of radiation detection. Specifically, this 
has motivated the search for new materials to improve the 
performance of room-temperature semiconductor radiation 
detectors.- - - Aluminum antimonide, which is a III-V semi- 
conductor with an indirect band gap of 1 .69 eV, is a promising 
candidates but its potential for radiation detection has not been 
conclusively established yet. 4,5 

The basic mechanism underlying radiation detection in a 
semiconductor can be summarized as follows: (1) an incom- 
ing quantum of radiation transfers energy to the electrons in 
the material creating electron-hole pairs, (2) the charge car- 
riers are separated in an externally applied electric field, (3) 
the magnitude of the resulting electric current is measured. 
To achieve maximum energy resolution, the "loss" of charge 
carriers due to scattering and recombination events must be 
minimized. This requires high carrier mobilities and very 
low concentrations of intrinsic defects and impurities. In ad- 
dition, a very high resistivity is desirable to maintain a low 
concentration of background free carriers and thereby to min- 
imize noise events. For instance, in germanium detectors that 
are used at cryogenic temperatures, defect and impurity con- 
centrations have to be reduced to less than 10 10 cm -3 and 
free carrier densities to less than 10 14 cm -3 to reach optimal 
performance. 6 To meet such stringent bounds a highly opti- 
mized manufacturing process is needed, which typically im- 
plies years of research and development. This situation ren- 
ders the screening of candidate materials extremely challeng- 
ing. Predictive computations, however, offer a possibility to 
circumvent this difficulty, since they allow studying the funda- 
mental limits of materials in a much more time effective and 
systematic fashion. Specifically, quantum mechanical meth- 
ods based on density-functional theory can be used to com- 



pute the formation energies of intrinsic and extrinsic defects 
and their equilibrium concentrations. They also allow one to 
obtain the defect-limited scattering rates, which determine the 
mobility and life time of free charge carriers. Since these cal- 
culations can be done for idealized systems, it is possible to 
separate the effects of different intrinsic and extrinsic defects 
and determine the ultimate performance limits. 

The objective of the present paper is to elucidate the role 
of extrinsic defects in aluminum antimonide on the basis of 
first-principles calculations. Specifically, we seek to estab- 
lish which elements are the most detrimental for the electronic 
properties of AlSb and which elements can be introduced in- 
tentionally to compensate for intrinsic defects and impurities. 
In short, it will be shown that (i) C, Si, and Ge substitute on 
Sb sites and act as acceptors with shallow transition levels, 
(ii) Sn can substitute on both Al and Sb sites and displays am- 
photeric characteristics, (in) S, Se, and Te substitute on Sb 
and act as donors with deep equilibrium transition levels, (iv) 
O is most likely to be incorporated as an interstitial and pre- 
dominantly acts as an acceptor, (v) P substitutes for Sb and is 
electrically inactive, (vi) C and O have the most detrimental 
effect on carrier transport, while (yii) Te is a weak scatterer 
and can therefore be used to manipulate the net concentra- 
tion of charge carriers with a modest effect on their mobili- 
ties. The present study represents a continuation of our earlier 
work on intrinsic defects in AlSbi and is connected to our re- 
cent fully quantum-mechanical calculations of defect-limited 
carrier mobilities. 8 

This paper is organized as follows. In Sec. JI] we give an 
overview of the methods and computational parameters. The 
formation energies of extrinsic defects and their binding en- 
ergies in complexes with intrinsic defects are summarized in 
Sec. [IllJ Then we discuss the effect of extrinsic defects on 
charge carrier mobility using a measure that is based on their 
associated lattice distortions in Sec. [iVl Finally, we analyze 
the impact of selected extrinsic defects on the net charge car- 
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FIG. 1 : (Color online) Overview of defect configurations considered in this study. Large violet, small gray, small orange and very small green 
spheres represent Sb, Al, impurity atoms, and H, respectively. 



rier and defect concentrations in the material and discuss op- 
timal doping schemes in Sec.lVl 



II. METHODOLOGY 



A. Defect thermodynamics 



ergy of a defect in charge state q can be written as 9 
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The equilibrium concentration, c, of a point defect depends 
on its free energy of formation, AG/, via 



c = cq exp (—AGf/hsT) . 



(1) 



where Co is the concentration of possible defect sites. The for- 
mation free energy AGf is usually approximated by the for- 
mation energy AEf which is legitimate if the vibrational en- 
tropy and the pressure- volume term are small^ In the present 
context we are interested in impurity-related defects in a bi- 
nary stoichiometric compound. In this case the formation en- 



where E^f is the total energy of the system containing the de- 
fect, rii denotes the number of atoms of type i, and /^ ulk is 
the chemical potential of component i in its reference state. 
Neglecting entropic contributions the chemical potentials of 
the reference phases can be replaced by the cohesive ener- 
gies at zero Kelvin. The formation energy depends on the 
chemical environment via A/iAisb which describes the vari- 
ation of the chemical potentials under different conditions. 
The range is constrained by the formation energy of AlSb, 
|A/iAisb| < AiJ/(AlSb), where for the present convention 
A/XAisb = -AiJ/(AlSb) and A/i A isb = Ai2/(AlSb) corre- 
spond to Al and Sb-rich conditions, respectively. The forma- 
tion energy also depends on the electron chemical potential, 
fjL e , which is measured with respect to the valence band max- 
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imum, i^vBM- The last sum in Eq. ([2]) runs over all the impu- 
rities in the system, where Ani denotes the difference in the 
number of atoms of impurity i between the ideal and the de- 
fective system, and A fa is the change in chemical potential 
of this component with respect to the bulk chemical poten- 
tial. For the latter we used the cohesive energy of the respec- 
tive most stable phase (e.g., O2 molecule, Si in the diamond 
structure, etc.) as computed within the present computational 
framework. 

The formation energies that are shown in Figs. H [6} and 
[7] have been obtained in the impurity-rich limit, i.e., A fa = 
eV. We have chosen this limit because it provides a good in- 
dication for how easy a certain impurity is incorporated during 
the growth process. In some cases, e.g., silicon and oxygen, 
this leads to negative formation energies. If one is interested 
in formation energies under different chemical conditions, the 
relevant formation energies can be easily obtained by shifting 
all formation upwards by the value of A fa chosen. In this 
way, one can e.g., determine the formation energies that are 
applicable if the AlSb sample is in contact with say a block of 
A1 2 3 ^ 

Given the formation energies for two different charge states, 
qi and q2, at the valence band maximum [fi e = in Eq. ©], 
we can also obtain the corresponding equilibrium (thermal) 
transition level, 

qi ~ Q2 

Note that equilibrium transition levels Ef r are distinct from 
electronic transition levels . While the former determine 
the equilibrium electron chemical potential at which a change 
in the charge state occurs, the latter are related to optical tran- 
sitions. In the present work, we obtained electronic transition 
levels from band structure calculations of defect cells. Exam- 
ples are shown in Fig. 01 

Binding energies for defect complexes are calculated ac- 
cording to 

AE b = A^ omplex - AE f , (4) 

where the sum runs over the individual defects from which the 
complex is constructed. Note that AE b changes with the elec- 
tron chemical potential, since AEf for both the complex and 
the individual defects depends on fa. According to Eq. © 
negative binding energies indicate a driving force for complex 
formation. 



B. Defect structures 

We consider C, Si, O, S, P, and H as potential impurities 
incorporated during the growth and include Ge, Sn, Se, and 
Te as possible compensating dopants. Figure [T] provides an 
overview of the defect configurations included in this work. 
To be specific, we took into account substitutional configura- 
tions with both tetrahedral (T^) and trigonal (Cs v ) symmetry 
as well as four variations of cation-cation bonded (CCBDX) 



configurations that have been discussed in connection with 
DX-centers in III-V and II- VI semiconductors (Addi- 
tional tests were carried out for several distorted substitutional 
configurations, including molecular dynamics simulations in 
order to check for the existence of other relevant local minima, 
but did not reveal any further configurations of importance). 
With regard to interstitial configurations, we included tetra- 
hedral (X^tet,Ai, ^Q,tet,Sb) and hexagonal sites (X^hex) as well 
as several split-interstitial configurations oriented along (100) 
and (110). The latter configurations were always observed to 
relax into a structure in which the bond between the impurity 
atom and its nearest neighbor is oriented slightly away from 
the (110) direction [compare Fig. |l(h)| . In the literature, this 
distorted structure has been previously referred to as "bridge" 
interstitial, and for simplicity in the following we adopt this 
convention. We typically calculated formation energies for 
charge states between q = — 2 and +2, and where necessary 
also included charge states —3 and +3. 

We also investigated the formation of defect complexes of 
group IV and VI elements with aluminum vacancies Vai and 
antimony antisites SbAi which are the two most important in- 
trinsic defects 1 



C. Computational details 

Calculations were performed within density-functional the- 
ory (DFT) using the local density approximation (LDA) and 
the projector augmented wave method 13,14 as implemented in 
the Vienna ab-initio simulation package^ - — 

In our earlier study 7 of intrinsic defects in AlSb we em- 
ployed a variety of supercell sizes in combination with a finite- 
size scaling procedure in order to mitigate the limitations of 
the supercell method. Given the large number of different ele- 
ments, defect structures, and charge states considered in the 
present work, such a comprehensive approach is, however, 
beyond our present computational capabilities. We therefore 
carried out the majority of the calculations for 64-atom super- 
cells, which based on the results of Ref. 7 implies that the error 
in our calculated formation energies due to the supercell ap- 
proximation is still only on the order of 0.1 eV. For Brillouin- 
zone integrations we used a 6 x 6 x 6 Monkhorst-Pack mesh. 19 
For a small set of configurations, we performed additional cal- 
culations using 216-atom cells (see Fig. [4]) and a 4 x 4 x 4 
Monkhorst-Pack fc-point mesh to verify the viability of our 
64-atom cell calculations. We used a plane- wave cutoff en- 
ergy of 500 eV, which is imposed by the pseudo-potential for 
oxygen, but for the sake of consistency was used for all cal- 
culations. A Gaussian smearing with a width of a = 0.1 eV 
was used to determine the occupation numbers. The ionic po- 
sitions were optimized until the maximum force was below 
20meV/A. For the charged defect calculations a homoge- 
neous background charge was added to ensure charge neutral- 
ity of the entire cell. 

For bulk AlSb we obtain a lattice constant of 6. 12 A in good 
agreement with the experimental room temperature value of 
6.13 A. The band gap is 1.12 eV, which is lower than the ex- 
perimental room temperature value of 1 .69 eV as expected for 
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a LDA calculation. A more detailed overview of the bulk 
properties of Al, Sb, and AlSb that are obtained using the 
present computational parameters can be found in Ref. 0. 

The underestimation of the band gap is a well-known de- 
ficiency of DFT within the LDA. In order to overcome this 
deficiency in Ref. [3 we have employed a post-mortem cor- 
rection scheme which relies on a separation of valence and 
conduction band states. Analysis of the electronic structure of 
the defects in this work reveals that in many cases the defect 
induces electronic levels in the band gap, which renders the 
separation into valence, conduction or localized states highly 
unreliable. Since as a result the application of the correction 
scheme employed in Ref. 7 becomes ambiguous, we have de- 
cided to abstain from any band gap corrections in the present 
work. Instead all results are presented and discussed with re- 
spect to the calculated band gap, as done previously by Du. 10 
For reference we have, however, included the experimental 
band gap in Figs. [2 O and[6] 

To compute the binding energies of defect complexes in- 
volving aluminum vacancies and antimony antisites, we used 
the defect formation energies obtained in Ref. 7 with for 64- 
atom cells and without band gap corrections 



D. Lattice distortion and effect on charge carrier mobilities 

Defects break the translational symmetry of an ideal crystal 
and therefore act as scattering centers for electrons and holes 
traveling through the crystal. The two most important contri- 
butions to the scattering cross section of a defect are related to 
the amount of charge localized at the defect (Coulomb scat- 
tering) and the lattice distortions that accompany defect relax- 
ation. Since we are primarily interested in comparing defects 
that have identical absolute charge states (e.g., C^ 1 vs Te^ 1 ), 
in the following we focus on the second contribution only. 
There are many different ways in which one can quantify lat- 
tice distortions. We have empirically found that a particularly 
useful measure is given by 

M 2 = (y dr\Vr(AV)\^J , (5) 

where AV denotes the potential difference AV between the 
defective and ideal cells. This quantity includes not only ionic 
relaxations but also changes in the electronic charge distribu- 
tion. 

In Eq. © we use the potential and not — for example — the 
charge density difference. This choice is motivated by a full 
fledged computation of the scattering rates 8 based on Fermi's 
golden rule for which one needs to evaluate matrix elements 
of the form 1 AV | **) j28i21 In that study^ we find that M 2 
reproduces very well the trends obtained from the scattering 
rate calculations based on Fermi's golden rule across a wide 
range of extrinsic defects. Since the computation of M 2 is, 
however, orders of magnitude more efficient, in the follow- 
ing we use the relative scattering strength M 2 to assess the 
impact of defects on carrier transport. We note that this analy- 
sis does not differentiate between electron and hole scattering. 



It is, however, unusual that a given defect has a significantly 
different impact on electron and hole scattering .- 

E. Self-consistent charge equilibration 

In Refs. and |22| we have described in detail a self- 
consistent scheme to determine the concentrations of defects 
and charge carriers as a function of temperature and chemical 
environment. Here, we only provide a brief summary of the 
most important equations. The concentrations of electrons, 
holes, and charged defects in a material are coupled to each 
other via the charge neutrality condition 

n e + n A = n h + n D , (6) 

where n e and rih are the electron and hole concentrations, re- 
spectively. They are given by 

POO 

n e = / D(E)f(E,ii e )dE (7) 

JCBM 

/VBM 
D(E)[l-f(E,n e )]dE, (8) 
-CO 

where D(E) denotes the density of states, f(E,fi e ) = 
{1 + exp [(E — n e ) /ksT}} -1 is the Fermi-Dirac distribu- 
tion, fi e is the electron chemical potential, VBM and CBM 
indicate the valence band maximum and the conduction band 
minimum, respectively. The concentrations of acceptors and 
donors, ua and ud, can be split into contributions due to in- 
trinsic and extrinsic defects 

n A = < + n e A xt and n D =rig + >nfg (9) 

The concentration of intrinsic acceptors (donors) is obtained 
by summing over all intrinsic defects with charge states q < 
(q > 0), 

acceptors donors 

i i 

where e is the unit of charge and the defect concentration is 
given by Eq. (Q]). 

The concentrations of extrinsic acceptors and donors are 
given by expressions that are analogous to Eq. (171)1) but sub- 
ject to the constraint that the total concentration of a given 
impurity or dopant is constant. A condition that can be easily 
achieved by proper normalization of the defect concentrations. 
The set of Eqs. (l6l-[TQl) can be solved self-consistently using an 
iterative algorithm. 

III. FORMATION ENERGIES AND DEFECT 
STRUCTURES 

A. Group IV elements: C, Si, Ge, and Sn 

The formation energies of defects involving the first four 
group IV elements (C, Si, Ge, and Sn) are shown in Fig. 0a- 
d). In all cases the most important configuration is the sub- 
stitutional defect on an Sb-site, Xsb (X=C, Si, Ge, or Sn). 
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FIG. 2: (Color online) Formation energies of impurity-related defects for (a-d) group IV (e) P, and (f-h) group VI elements calculated as a 
function of the electron chemical potential with A/ii mp = OeV in Eq. (|2]). Symbols indicate equilibrium transition levels [see Eq. (|3)]. The 
slope of the lines corresponds to the charge state. The white region indicates the calculated band gap while the lighter gray region indicates 
the difference between the experimental and the calculated band gap. 



Regardless of the initial condition, we find the defects to re- 
lax into a configuration with Td symmetry. The Xsb defects 
occur predominantly in charge state q = — 1 and thus act as 
single acceptors. This observation is in agreement with Hall 
and resistivity measurements.— - — 

Within the error bars of our calculations the 0/ — 1 equi- 
librium transition levels for C, Si, and Ge coincide with the 
VBM whereas for Sn the equilibrium transition level is lo- 
cated 0.1 eV above the VBM (compare Fig.O. In agreement 
with our calculations, shallow acceptor levels have been ob- 



served for Si by Bennett et al. while C has been found to 
act an acceptor by McCluskey et al. & For Xsb defects, our 
calculations do not reveal any electronic states inside the band 
gap that could give rise to strong optical signatures. 

For C the substitutional defect on an Al-site Xau which 
always adopts a local Td symmetry, is very high in energy 
and therefore should not occur in thermodynamic equilibrium. 
For Si, Ge, and Sn, however, this defect becomes increas- 
ingly more important. In fact in the case of Sn, the forma- 
tion energies for substitutional defects on Al and Sb-sites are 
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FIG. 3: (Color online) Equilibrium (left) and electronic (right) transition levels of the most relevant impurity related defects for group IV and 
VI elements. The equilibrium transition levels have been extracted from the formation energies using Eq. (|3j and correspond to the position of 
the symbols in Figs.[2]and[5] The electronic transition levels have been identified using band structure calculations of the type shown in Fig. [4] 
0sb,i and 0sb,2 stand for 0sb(C3v) and 0sb(/3 — CCBDX), respectively. All other substitutional defects have Td symmetry. The oxygen 
interstitial has tetrahedral symmetry. 
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FIG. 4: (Color online) Band diagrams for (a) an ideal system as 
well as for systems containing (b) a GeAi defect and (c) a 0sb (/?- 
CCBDX) defect. Conduction and valence bands are indicated by red 
and blue lines, respectively. Defect states are shown in thick yellow 
lines. All calculations were carried for 216-atom supercells in charge 
state q = 0. 



configuration for Si, Ge, and Sn. All interstitials are, how- 
ever, significantly higher in energy than the most stable sub- 
stitutional defects whence in thermodynamic equilibrium their 
concentrations should be negligible. Since the migration bar- 
riers for substitutional defects are, however, typically very 
large, it is possible that interstitial configurations play a role 
with respect to the diffusion of group IV elements. 

We also considered defect complexes of group IV impu- 
rities with the intrinsic defects Vai and SbAi- The formation 
energies for Xsb — Vai defects which act as acceptors are 
rather high due to the like charges of the two isolated de- 
fects. The Xsb — SbAi complexes on the other hand display 
rather small formation energies that are typically less than 
1 eV higher than the formation energy of the most stable con- 
figuration. These defects exhibit amphoteric behavior with the 
donor state dominating over the wider range of the band gap 
[also see Fig. Ha)]. 



very similar. These results are consistent with the experiments 
by Shaw2£ who observed that Sn can occupy both Al and Sb 
sites. For Si, Ge, and Sn the calculations predict Xai to be a 
donor with a deep +1/0 equilibrium transition level (compare 
Fig. [3]). For Si, Ge, and Sn we also find electronic levels lo- 
cated inside the band gap. This is explicitly demonstrated for 
the case of GeAi in Fig.[4tb), which reveals a localized defect 
state approximately 0.24 eV below the CBM. 

The similar formation energies of SnAi donors and Snsb ac- 
ceptors gives rise to an interesting type of amphoteric behavior 
for Sn-doped AlSb. Usually this property is associated with a 
single defect type, which undergoes electronic (and usually 
also structural) changes as a function of charge state. In the 
case of AlSb:Sn, however, the amphoteric behavior is due to 
two different defects occupying two different lattice sites that 
possess opposite electronic characteristics. This leads to an 
intriguing opportunity for compensation doping that will be 
explored in Sec. lVDl 

Regarding interstitial defects, both tetrahedral configura- 
tions X^ tet) Ai and X^ tet ,sb [Figs. [T(c)| and |l(d)| are found to 
be local minima on the energy landscape with Xi t e t ai 

be- 
ing slightly more stable. The hexagonal interstitial X^hex is 
a local minimum for C only, whereas it relaxes to the X^ tet ,Ai 



B. Group V elements: phosphorus 

Phosphorus is included in the present work as a prototypi- 
cal example for group V elements which are isoelectronic to 
Sb. The calculated formation energies are shown in Fig.[2e) 
which shows that P prefers substitution on a Sb-site maintain- 
ing Td symmetry. Psb is neutral throughout the band gap and 
thus does not introduce any transition levels inside the gap. 
All other possible configurations are much higher in energy. 

C. Group VI elements: S, Se, and Te 

Since oxygen behaves very different from the other group 
VI elements we first consider the less complex cases of S, Se, 
and Te. A detailed discussion of oxygen is given in the next 
section. 

For S, Se, and Te the most stable defect configuration for 
electron chemical potentials over the widest range of the band 
gap is Xsb which possesses Td symmetry and acts as a donor 
[Figs. Etf-h)]- The +1/0 equilibrium transition levels for 
these defect are deep for all three elements (Fig. [3}. While all 
interstitial configurations considered are local minima, they 
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FIG. 5: (Color online) Formation energies of oxygen-related defects under Al-rich and Sb-rich conditions [A//0 
indicate equilibrium transition levels. The slope of the lines corresponds to the charge state. 
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are always much higher in energy than the respective Xsb de- 
fect. Although there is supporting experimental evidence of 
the existence of shallow donor levels r&2Z it is firmly estab- 
lished that the ground- state of these donors for electron chem- 
ical potentials near the conduction band minimum is actually 
a DX-center^&2? associated with a deep level^"— As indi- 
cated in Sec. Ill Bl we have investigated a variety of possible 
DX-center structures. While we find several of these config- 
urations to be local energy minima for charge state q = — 1, 
their formation energies are always higher than the respec- 
tive Td configurations in charge state as long as the electron 
chemical potential remains within the band gap. This behav- 
ior can be attributed to the LDA which in this instance fails 
to properly capture the energetics of localized vs delocalized 
defect states. A detailed discussion is, however, beyond the 
scope of this article and will be published elsewhere. 

Defect complexes of X$b with SbAi antisites for the group 
VI impurities show similar trends as in the case of the group 
IV impurities being about 1 eV higher than the isolated X$b 
defects [see Fig.[7tc)]. In all cases, these complexes display 
donor character. In contrast to the group IV elements, how- 
ever, the Xsb — Vai complexes have much lower formation 
energies resulting from the Coulombic interaction of the op- 
positely charged isolated defects. In fact, for electron chem- 
ical potentials in the upper half of the band gap, this mutual 
attraction renders the Xsb — Vai acceptor complex the most 
stable defect. The donor character of X$b combined with the 
acceptor character of Xsb — Vai could allow for an interest- 
ing self-compensation effect. Since this would require that 
aluminum vacancies remain mobile down to very low temper- 
atures (< 400 K) in order to achieve full equilibrium, such a 
scenario is, however, unlikely to occur in reality. 



D. Group VI elements: oxygen 

Oxygen behaves distinctly different from all the other ele- 
ments considered so far , in that its interstitial configurations 
have very low formation energies and there is a multitude of 
different defect configurations which are all very close in en- 
ergy (Fig.0. 



Over the widest range of the electron chemical potential, 
the most stable defect is the tetrahedral interstitial surrounded 
by Al atoms 0i, t et,Ai in charge state q = — 2 [Fig. 1 1(c)) . This 
acceptor defect displays transition levels 0.1 eV (0/ — 1) and 
0.5 eV (-1/ - 2) above the VBM (see Fig.0. 

For Al-rich conditions and electron chemical potentials in 
the lower half of the band gap, the lowest energy defect is 
an interstitial which assumes a bridgelike structure [Fig. 1 1(h)) 
and is indicated as 0^bb in Fig. [5] It maintains the neu- 
tral charge state throughout the band gap. Very close to the 
VBM the &sb(C3 V ) defect emerges as the ground-state con- 
figuration. This defect adopts a configuration with C$ v sym- 
metry and is associated with a pronounced lattice distortion 
[Fig. |l(a)| . It displays amphoteric characteristics with equi- 
librium transition levels 0.2 eV (+1/0) and 0.4 eV (0/ - 1) 
above the VBM. As shown in Figs. [3] and [He), it also intro- 
duces two electronic levels into the band gap that are located 
0.45 eV above the VBM and right at the CBM, respectively. 

For Sb-rich conditions and electron chemical potentials in 
the lower half of the band gap the 0sb — SbAi defect complex 
which exhibits donor character prevail. 

Finally, there are several other configurations with forma- 
tion energies within about 0.3 eV of the respective ground- 
state configurations. These include the (0 — Sb)^ 00 ^ split- 
interstitial [Fig. |l(g)| , which is neutral throughout the band 
gap and electrically inactive, as well as the 0sb(/5 — CCBDX) 
defect [Fig. [T(m)) . 

We observe a large number of oxygen configurations with 
very similar formation energies, and most of these defects 
lead to one or more deep equilibrium transition levels. Deep 
states induce defect-mediated recombination and/or trapping 
and thereby limit the lifetime of charge carriers. Oxygen in- 
corporation is therefore very detrimental to the operation and 
energy resolution of a detector device. 



E. Hydrogen 

For hydrogen, we focused on interstitial configurations and 
considered both tetrahedral sites and the bridge bond config- 
uration [Fig. |l(i)| . The results are shown in Fig. [6) For elec- 
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FIG. 6: (Color online) Formation energies of hydrogen-related de- 
fects (A /in = OeV in Eq. ©). Symbols indicate equilibrium tran- 
sition levels. The slope of the lines corresponds to the charge state 
[see Eq. ©]. 



tron chemical potentials in the lower half of the band gap, the 
bridge bond configuration is the energetically most favorable 
one. It acts as a deep donor with an equilibrium transition 
level 0.3 eV below the CBM. For electron chemical potentials 
in the upper half of the band gap, the H^ tet) Ai defect prevails 
which behaves as an acceptor with an equilibrium transition 
level 0.1 eV above the VBM. We also note that although the 
formation energies for H^bb and H^ tet ,Ai in the neutral charge 
state are virtually identical, the structures are nonetheless dis- 
tinct. 



F. Defect complexes 

In this section, we consider the interaction of the two most 
important intrinsic defects 7 SbAi and Vai with group IV and 
VI impurities. In the preceding sections, we have already seen 
that these defect complexes can have very low formation ener- 
gies, in particular in the case of group VI elements, where for 
certain electron chemical potentials they are the most stable 
defect configuration. 

In Fig. [71 we show the binding energies and the respective 
formation energies of Sn and Te-complexes that are prototyp- 
ical for group IV and VI elements, respectively. To illustrate 
the trends we show in Fig. [8] the most negative binding ener- 
gies (indicating the strongest driving force for defect complex 
formation) for all group IV and VI impurities considered. 

In the case of Sn, both types of complexes display strong 
binding for electron chemical potentials in the lower two- 
thirds of the band gap with binding energies as negative as 
— 2.2 eV for Snsb — SbAi at the valence band maximum. On the 
other hand, for electron chemical potentials in the upper third 
of the band gap, the interaction is strongly repulsive. These 
trends hold for the other group IV elements as well with the 
exception of the Csb — Vai complex which is never bound. 
In general, the vacancy complexes have high formation ener- 
gies and act as acceptors, while the antisite complexes have 
rather low formation energies near the valence band edge and 
display mostly donor characteristics. 

As exemplified by Te, the trends described above for X$b — 



SbAi prevails also for group VI elements. With respect to 
-^sb — Va\ complexes with group VI impurities, we obtain 
binding energies which are consistently negative throughout 
the band gap. As discussed in Sec. MI CI the vacancy com- 
plexes act as acceptors, have relatively low formation ener- 
gies, and are actually the most stable defect for electron chem- 
ical potentials near the CBM. 



IV. RELATIVE SCATTERING STRENGTHS 

The relative carrier scattering strengths calculated using 
Eq. ([5]) for the most important intrinsic and extrinsic defects 
vary over several orders of magnitude as shown in Fig. [9j 
While the absolute values of the relative scattering strength 
bear no meaning, they still allow us to compare different ex- 
trinsic and intrinsic defects with respect to their impact on 
charge carrier mobility. 8 

Figure [9] clearly shows that oxygen induces by far the 
strongest carrier scattering and thus has the most detrimental 
impact on charge carrier mobilities, with carbon also show- 
ing a very high scattering strength. The relative scattering 
strengths for these impurities are about one order of magni- 
tude larger than for the aluminum vacancy, which is the most 
detrimental intrinsic defect. 

For the elements within the same group of the periodic ta- 
ble, the scattering strength generally decreases with increas- 
ing atomic number. This effect is owed to the diminishing 
size mismatch between the impurity atom and Sb, for which 
most impurities substitute, leading to decreased lattice distor- 
tion upon substitution. As a result, the elements that have 
the least negative effect on carrier mobility, are Sn, Se, and 
particularly Te. In fact, many of the substitutional impurities 
scatter carriers less strongly than the intrinsic defects. These 
elements are therefore excellent candidates for doping the ma- 
terial in order to achieve carrier compensation with minimal 
impact on mobility, an opportunity that is explored in the fol- 
lowing section. 



V. CHARGE CARRIER CONCENTRATIONS 

In the previous sections we have determined the formation 
energies of a wide range of extrinsic defects and discussed 
which of them act as donors or acceptors. Furthermore, we 
have obtained their relative scattering rates, which enabled us 
to identify the defects which are, respectively, the least and 
most detrimental to the electron and hole mobilities in the 
material. Now, we are in a position to combine these data 
to predict net charge carrier concentrations and to determine 
conditions for obtaining optimal material. For radiation de- 
tection, the objective is to obtain a material which features 
maximal charge carrier mobilities, which requires small scat- 
tering rates, and very high resistivity, which requires low net 
charge carrier concentrations. 

In the analysis that follows, we consider four important 
cases. First, we analyze the pure material (containing only in- 
trinsic defects). Then, we include in our model unintentional 
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FIG. 7: (Color online) (a,c) Formation and (b,d) binding energies of defect complexes of intrinsic defects with (a,b) Sn and (c,d) Te. The two 
cases are exemplifying the situation for group IV and VI elements, respectively. 
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FIG. 8: (Color online) Most negative binding energies for Xsb — Vai 
and Xsb — SbAi complexes where X is a group IV or VI element. 
Following Eq. $4$ negative binding energies indicate a driving force 
for defect complex formation. 



charge equilibration scheme to study the behavior of extrinsi- 
cally doped AlSb, we therefore have to discuss to which ex- 
tent the band gap error influences the results. According to 
Eq. (|7]) the intrinsic electron and hole concentrations depend 
exponentially on the band gap. To be consistent with the cal- 
culations discussed above, we now have to use the LDA band 
gap and the corresponding formation energies for the intrin- 
sic defects. The absolute carrier concentrations are therefore 
no longer comparable to our calculations in Ref . [J For the 
present purpose, we are, however, interested not in absolute 
carrier concentrations but in net charge carrier concentrations, 
i.e. the difference between the concentrations of electron and 
holes, 



An : 



(11) 



impurities from the growth process, followed by analysis of 
carrier compensation with intentional Te doping. Finally, we 
consider the case of Sn doping, which exhibits unique behav- 
ior due to its amphoteric character. 



A. Pure material 

In Ref. 7, we employed the self-consistent charge equilibra- 
tion scheme described in Sec. lIIEI to determine the concentra- 
tions of intrinsic defects and charge carriers in pure AlSb (no 
impurities) using defect formation energies calculated from 
DFT, which had been corrected for the band gap error. As 
discussed in Sec. Ill CI in the present work we did not employ 
any such corrections. Before we can apply the self-consistent 



The impact of the band gap error on An is not immediately 
obvious. In Fig. [I0l we therefore compare the net electron 
concentrations for pure AlSb using as-calculated formation 
energies as well as formation energies, which have been cor- 
rected for the band gap error. We find that while the actual val- 
ues deviate somewhat, the conduction type (n vs p), chemical 
potential dependence, and the order of magnitude of An are 
very similar. This finding enables us to use the as-calculated 
formation energies to compute net charge carrier concentra- 
tions for AlSb with extrinsic defects. Figure [TO] also shows 
that pure AlSb exhibits n-type behavior over the entire range 
of the chemical potential difference. In contrast, as-grown 
AlSb has been found to exhibit p-type conductivity. In the 
following section, we show that this apparent discrepancy is 
caused by the accidental incorporation of impurities during 
growth. 
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FIG. 10: Net electron concentration for pure AlSb at three different 
temperatures using band gap corrected (solid lines) and as-calculated 
(dashed lines) formation energies assuming full equilibrium. The 
horizontal bar indicates the range of chemical potentials that are ac- 
cessible based on the calculated compound formation energy. 



B. Unintentional background doping 

In the experimental growth process, certain impurities are 
impossible to avoid. In particular this includes oxygen and 
carbon Oxygen can be incorporated from the gas phase 
or the components of the growth chamber and is also present 
in the raw starting materials. Carbon enters because graphite 
susceptors are used to heat the melt. 35 Using a highly opti- 
mized growth process, it is possible to reduce the concentra- 
tion of oxygen to about 10 15 cm -3 and the concentration of 
carbon to about 10 17 cm~ 3 ^ 

In Fig. [TT] we show the net charge carrier concentrations 



cm 



for AlSb with c(C) = 10 17 cm" 3 and c(0) = 10 
In fact, these concentrations cause the conductivity to switch 
from n-type to p-type, which resolves the apparent discrep- 
ancy between theory and experiment mentioned above. The 
inversion of the conductivity-type is indeed expected since in 
SecJIIlboth of these impurities have been identified as accep- 



tors. 



C. Compensation with Te 

Since the background doping by C and O renders the ma- 
terial p-type, compensation requires addition of donor ions. 
To simultaneously achieve high carrier mobilities and low net 
carrier concentrations, the dopant should possess a scattering 
potential which is as small as possible. According to Figs. [2] 
and [9] Te is an ideal candidate based on these criteria. In 
Fig. [TT] we show the effect of adding Te at a concentration 
of 10 17 cm -3 to AlSb containing C and O impurities. We find 
that by using Te doping in this concentration range the net 
charge carrier concentration can be tuned very efficiently. The 
net carrier concentration changes as a function of the chemical 
potential difference between Al and Sb, with a minimum close 
to the Sb-rich limit (typical annealing condition). Intentional 
doping with Te therefore enables us to compensate the back- 
ground doping due to C and O impurities while minimizing 
the detrimental impact on the carrier mobilities. In fact, using 
these data as a guidance for the experimental growth process 
led to a dramatic improvement of the material for radiation 
detection, which is reported elsewhere.— 



D. Ambipolar doping with Sn 

According to Fig. [2 Sn can be incorporated both on the Sb 
and the Al sublattice, which gives rise to acceptor and donor- 
type defects, respectively. Also, Fig. [9] shows that both of 
these defects are weak scatterers compared to the aluminum 
vacancy. Figure 02 a) illustrates the effect of Sn on the net 
charge carrier concentration. In thermodynamic equilibrium 
and for a Sn concentration of about 10 17 cm -3 , the material is 
n-type if grown under Sb-rich conditions and p-type if grown 
under Al-rich conditions [case II in Fig. H2a)]. In the pres- 
ence of C and O impurities this behavior is retained although 
the transition point is shifted toward Sb-rich conditions [case 
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FIG. 11: Net charge carrier concentration for (I) pure AlSb, (II) 
AlSbwithc(C) = 10 17 cm- 3 and c(0) = 10 1 
AlSb containing c(C) = 10 17 cm" 3 , c(0) = 
c(Te) = 10 17 cm -3 for two different temperatures assuming full 
equilibration and using as-calculated formation energies. If An in 
Eq. dTTb is positive (negative) as shown by the solid (dashed) lines the 
material exhibits n-type (p-type) conductivity. The horizontal bar in- 
dicates the range of chemical potentials that are accessible based on 
the calculated compound formation energy. 



C in Fig. [Eta)]. By increasing the Sn concentration the mini- 
mum in An can be shifted back toward the impurity-free case 
[case IV in Fig. [T2t aYl. Figure \Y2\ b) reveals that the origin 
of this behavior is the subtle equilibrium between Sn^ 1 and 
Sn^ 1 defects, which act as acceptors and donors, respectively. 

We suggest that this intriguing property of Sn-doped AlSb 
can be exploited to achieve self-compensation. In a typical 
experiment, an AlSb single crystal is grown under Al-rich 
conditions close to the melting point of 1330 K. The growth 
process is sufficiently slow to assume that we are near ther- 
modynamic equilibrium, which implies that the vast majority 
of Sn occupies Sb- sites. The growth step is usually followed 
by a high-temperature anneal at about 1200K under Sb-rich 
conditions. This inversion of chemical conditions leads to a 
strong driving force for redistributing Sn from the Sb to the Al 
sublattice. This conversion can be mediated by aluminum va- 
cancies which at this temperature are both very mobile^ and 
exist in high concentrations. 7 By carefully controlling the an- 
nealing time it should thus be possible to achieve a Snsb / SnAi 
ratio close to one. At this point the material would be perfectly 
compensated due to the opposing effects of Snsb and SnAi de- 
fects. Due to the moderate scattering strength of Sn-related 
defects such a material would exhibit high carrier mobilities 
in conjunction with high resistivity, as needed for the radiation 
detection applications. 
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FIG. 12: (a) Net charge carrier concentration for (I) pure AlSb, (II) 



AlSb with c(Sn) = 10 17 cm" 3 , (III) AlSb with c(C) 
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with c(C) = 10 17 cm- 3 , c(0) = 10 15 cm- 3 , and c(Sn) = 2 x 
10 17 cm -3 at a temperature of 1200 K assuming full equilibration 
and using as-calculated formation energies. Note that in cases (II- 
IV) going from Sb-rich to Al-rich conditions changes the conductiv- 
ity type from n top. (b) Defect concentrations for case (II) in (a) 
revealing that the inversion of the conductivity type observed in Sn- 
doped AlSb is caused by a shift from Sn^j 1 to Sn^ 1 . The horizontal 
bars indicates the range of chemical potentials which are accessible 
based on the calculated compound formation energy. 



VI. CONCLUSIONS 

In this work we have investigated the thermodynamic and 
electronic properties of group IV and group VI impurities as 
well as P and H in AlSb using density-functional theory cal- 
culations. First, we computed the formation energies for a 
wide range of possible defect configurations and charge states 
to identify the most important defect configurations, their 
electronic characteristics, and their equilibrium and electronic 
transition levels. We found that C, Si, and Ge preferentially 
substitute for Sb, act as acceptors and exhibit shallow tran- 
sition levels. Sn, however, can substitute for Sb and Al and 
displays amphoteric behavior. 

With regard to group VI elements, S, Se, and Te substitute 
for Sb and have deep donor levels. Oxygen in contrast can be 
incorporated in several different configurations that display a 
variety of electronic characteristics and transition levels, pre- 
dominantly of the acceptor type. For the widest range of elec- 
tron chemical potentials a tetrahedral interstitial configuration 
is the lowest in energy. 
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Phosphorus substitutes for Sb and is electrically inactive 
while H is incorporated on interstitial sites and can act as ac- 
ceptor or donor depending on conditions. 

With respect to defect complexes, our calculations show 
that group IV elements bind to both aluminum vacancies (Vai) 
and Sb antisites (SbAi) for electron chemical potentials in the 
lower two-thirds of the band gap. (The Csb — Vai complex is an 
exception and is never bound). The vacancy complexes have 
high formation energies and act as acceptors while the anti- 
site complexes have rather low formation energies near the 
valence band edge and display mostly donor characteristics. 
Similar trends hold for group VI elements but in this case the 
vacancy complexes show negative binding energies over the 
entire range of the band gap. In addition, the formation en- 
ergies of antisite complexes are very low near the conduction 
band edge, which for S, Se, and Te renders them even more 
stable than the individual Xsb defects in a small range of elec- 
tron chemical potentials in the upper third of the band gap. 

To characterize the effect of different impurities on the 
charge carrier mobilities, we calculated relative scattering 
strengths. This analysis revealed that C and O have by far 
the most detrimental impact on carrier transport. The scat- 
tering strengths for Sn, Se, and Te, however, are similar or 
even smaller than for the most important intrinsic defects. Tel- 
lurium, in particular, is therefore an excellent candidate for 
compensation doping the material. 



The interplay of electronic characteristics and impurity 
concentrations was investigated using a self-consistent charge 
equilibration scheme that allowed us to compute net charge 
carrier concentrations. While pure material was found to ex- 
hibit n-type conductivity, the incorporation of C and O leads 
to p-type material. This observation is consistent with ex- 
perimental measurements which show p-type conductivity in 
as-grown samples which contain an unintentional background 
concentration of C and O impurities. Our calculations fur- 
thermore revealed how Te can be used to compensate these 
impurities in order to achieve minimal net charge carrier con- 
centrations. Finally, we discussed how the amphoteric charac- 
ter of Sn could be used to achieve self-compensated material 
via controlled annealing. 
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